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Abstract 


Variations in fi, the total density of the Universe, leave a clear and distinctive 
imprint on the power spectrum of temperature fluctuations in the cosmic mi¬ 
crowave background (CMB). This signature is virtually independent of other 
cosmological parameters or details of particular cosmological models. We 
evaluate the precision with which 17 can be determined by a CMB map as a 
function of sky coverage, pixel noise, and beam size. For example, assuming 
only that the primordial density perturbations were adiabatic and with no 
prior information on the values of any other cosmological parameters, a full- 
sky CMB map at 0.5° angular resolution and a noise level of 15 /uK per pixel 
can determine 17 with a variance of 5%. If all other cosmological parameters 
are fixed, 17 can be measured to better than 1%. 
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Determination of the geometry of the Universe remains perhaps the most compelling 
problem in cosmology. Alternatively stated, what is the mean total energy density of the 
Universe? The answer to this question will reveal the ultimate fate of the Universe. If the 
density D (in units of the critical density Pc = SHq/SttG, where Hq is the Hubble constant) 
is greater than unity, the Universe is closed and will eventually recollapse; if it is less than 
unity, the Universe will expand forever; and if D = 1, the expansion will asymptotically 
decelerate to zero. 

Theoretical considerations favor a critical (D = 1) Universe, and inflation provides a 
generic mechanism for obtaining D = 1. However, luminous matter provides less than one 
percent of this mass. Various inferences of D by dynamical means have hinted at substantial 
amounts of unseen mass, but most traditional methods of determining D are plagued by 
systematic uncertainties. Furthermore, virtually all dynamical methods of obtaining D give 
the mean density in only nonrelativistic matter, and thus cannot discriminate between an 
open Universe and a flat Universe that is dominated by vacuum energy (i.e., a cosmological 
constant). 

Recently, it was proposed that temperature anisotropies in the cosmic microwave back¬ 
ground (CMB) might be used to determine the geometry of the Universe [Q. Features 
(known as “Doppler peaks,” or more accurately as acoustic peaks) in the CMB angular 
power spectrum result from acoustic oscillations in the photon-baryon fluid before the pho¬ 
tons decouple. The characteristic wavelength of these fluctuations is the sound horizon at 
decoupling (the distance an acoustic disturbance propagates from t = 0 until decoupling), 
which subtends an angular scale on the sky today of 0 ~ The dependence on D arises 

directly from the geometry of the Universe, and this angular scale is largely independent 
of other cosmological parameters. Thus, the location of the hrst Doppler peak provides a 
robust determination of D. A CMB map with hne angular resolution also constrains the 
other cosmological parameters by measuring the angular locations and amplitudes of the 
higher Doppler peaks. 

In this paper, we evaluate the precision with which D can be determined with high- 
resolution CMB maps [^. We work within the context of models with adiabatic primordial 
density perturbations, although similar arguments apply to isocurvature models as well [Q, 
and we expect the power spectrum to distinguish clearly the two classes of models. We also 
briefly consider what information on other cosmological parameters the CMB can provide. 

A given cosmological theory makes a statistical prediction about the distribution of CMB 
temperature fluctuations, expressed by the angular power spectrum 




/AT(m) AT(n) 
\~ro 7^ 


^ mn=cos 9 


J2^^CePe{cos9), 


( 1 ) 


where AT(h)/To is the fractional temperature fluctuation in the direction h. Pi are the 
Legendre polynomials, and the brackets represent an ensemble average over all observers 
and directions. The mean CMB temperature is Tq = 2.726 ± 0.010 K IQ. Since we can only 
observe a single microwave sky, the observed multipole moments will be distributed 
about the mean value Ci with a “cosmic variance” oi ~ ^2/(2^ -f- l)Cp no measurement 
can determine the Ci to better accuracy than this variance. 

We consider an experiment which maps a fraction /sky of the sky with a gaussian beam 
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with full width at half maximum 6*f„hm and a pixel noise Upix = where s is the de¬ 

tector sensitivity and tpix is the time spent observing each 6*f„hm x ^'fwhm pixel. We adopt the 
inverse weight per solid angle, w~^ = (crpix6'fwhm/7o)^, as a measure of noise that is pixel-size 
independent p . Current state-of-the-art detectors achieve sensitivities of s = 200 /iK -y/sec, 
corresponding to an inverse weight of w~^ ~ 2 x 10“^® for a one-year experiment. Re¬ 
alistically, however, foregrounds and other systematic effects may increase the noise level; 
conservatively, w~^ will likely fall in the range (0.9 — 4) x 10“^“^. Treating the pixel noise as 
gaussian and ignoring any correlations between pixels, estimates of Ci can be approximated 
as normal distributions with a variance (modihed from Ref. ||^) 
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Given a spectrum of primordial density perturbations, the Ci are obtained by solving 
the coupled equations for the evolution of perturbations to the spacetime metric and per¬ 
turbations to the phase-space densities of all particle species in the Universe. We consider 
models with initial adiabatic density perturbations hlled with photons, neutrinos, baryons, 
and collisionless dark matter; this includes all inflation-based models. We begin with ap¬ 
proximate analytic solutions for the scalar and tensor metric perturbations. Our 
calculation includes polarization |^, scale dependence of the initial perturbation spectrum 
1^, and the large-angle integrated Sachs-Wolfe effect from a cosmological constant m- To 
a good approximation, reionization can be parameterized by the optical depth r to the sur¬ 
face of last scatter |l|; anisotropies on scales much smaller than the horizon at reionization 
are suppressed by while those on larger scales are unaffected. The geometry of the 
Universe is then accounted for by shifting the moments, = 1) |l|, and 

approximating the large-angle integrated Sachs-Wolfe effect in an open universe |]^ . We do 
not here account for massive neutrinos (hot- or mixed- dark-matter models), but the power 
spectrum is altered only slightly by trading some of the nonrelativistic matter for neutrinos 
and our results should be unchanged 0. 

The CMB power spectrum depends upon many parameters. In the present anal¬ 
ysis we include the following set: the total density hi; the Hubble constant, Hq = 
100 hkmsec”^ Mpc“^; the density of baryons in units of the critical density, the 

cosmological constant in units of the critical density. A; the power-law indices of the initial 
scalar- and tensor-perturbation spectra, ns and n^; the amplitudes of the scalar and tensor 
spectra, parameterized by Q, the total CMB quadrupole moment, and r = Qt/Qs, the 
ratio of the tensor and scalar quadrupole moments; the optical depth to the surface of last 
scatter, r; the deviation from scale invariance of the scalar perturbations, Orun = dn/dlnk] 
and the effective number of light-neutrino species at decoupling, W- Thus for any given set 
of cosmological parameters s = {H, Qhh'^, h, ns, A, r, n^, drum t, Q, we can calculate the 
mean multipole moments ^^(s). 

We now wish to determine the capability of CMB maps to determine these cosmological 
parameters. The answer to this question will depend on the measurement errors cx/, and on 
the underlying cosmological theory. If the actual parameters describing the Universe are Sq, 
then the probability distribution for observing a CMB power spectrum which is best £t by 
the parameters s is 
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-(s-So) • H • (s-So) 


F(s) oc exp 

where the curvature matrix [a] is given approximately by 




dCi{so) 9Q(so) 


dsi 
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with a£ as given in Eq. The covariance matrix [C] = [a]“^ is an estimate of the stan¬ 
dard errors that would be obtained from a maximum-likelihood £t to data: the variance 
in measuring the parameter Si (obtained by integrating over all the other parameters) is 
approximately CJ . If some of the parameters are known, then the covariance matrix for 
the others is determined by inverting the submatrix of the undetermined parameters. For 
example, if all parameters are hxed except for Sj, the variance in Sj is simply In 

previous work, variances were estimated for small subsets of the parameters with Monte 
Carlo calculations 


131; the present approach can be used to reproduce these results. 


Fig. I displays the variance in as a function of the beam width 0fwhm for different 
noise levels and for /gky = 1. For different values of /sky, replace w —>• wfsky and scale by 
/sky'^^ [c.f., Eq. (H)]. The underlying model assumed here for the purpose of illustration is 
“standard CDM,” given by s = {1, 0.01, 0.5,1, 0, 0, 0, 0, 0, Qcobe, 3}, where Qcobe = 20 fiK 
is the COBE normalization |^. The solid curves show the Cqq obtained by inversion of 
the full 11 X 11 curvature matrix [a] for w~^ = 2 x 10“^^, 9 x 10“^^, and 4 x lO”^"^. These 
are the sensitivities that can be attained at the given noise levels with the assumption of 
uniform priors (that is, including no information about any parameter values from other 
observations). The dotted curves show the obtained by inversion of the Vt-Q submatrix 
of [a]; this is the variance in that could be obtained if all other parameters except the 
normalization were hxed, either from other observations or by assumption. Realistically, the 
precision obtained will fall somewhere between these two sets of curves. Other underlying 
models, including low-0 models, give similar sensitivities. Although parameters other than 
O will have some weak effect on the position of the hrst Doppler peak, they will also alter 
the power spectrum at smaller angular scales. Therefore, the higher multipole moments 
accessible with smaller beam widths help constrain the other parameters and make the 
determination of O from the location of the hrst Doppler peak more precise. 

Early reionization tends to wash out the structure of the power-spectrum features, de¬ 
creasing the precision of the parameter estimates. To illustrate this ehect, the curves in 
Fig. I show the same results as in Fig. |I|, but for a reionized model in which r = 0.5. As 
expected, the sensitivity to D decreases, although it remains signihcant even for r as large 
as one half. 

At £ ^ 1000, nonlinear ehects become signihcant and linear power-spectrum calculations 
become unreliable. Therefore, we extend the sum in Eq. (^ only up to £ = 1000. With 
improved calculations, the sensitivities at small beam widths could conceivably be improved. 

We have also investigated the sensitivity of CMB mapping experiments to the other 
cosmological parameters listed above. Our results suggest that a map with 0.5° angular 
resolution may also provide interesting constraints to A with minimal assumptions, and to 
the other parameters with reasonable priors. In particular, the experiments should be able 
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FIG. 1. The variance on H that can be obtained with a full-sky mapping experiment as a 
function of the beam width 6*fw hm for noise levels = 2 x 10“^®, 9 x 10“^®, and 4 x 10“^^ 
(from lower to upper curves). The underlying model is “standard CDM.” The solid curves are the 
sensitivities attainable with no prior assumptions about the values of any of the other cosmological 
parameters. The dotted curves are the sensitivities that would be attainable assuming that all 
other cosmological parameters, except the normalization, were hxed. The results for a mapping 
experiment which covers only a fraction /sky of the sky can be obtained by replacing w —> tc/sky 
and scaling by [c.f., Eq. (2)]. 


to distinguish between a flat matter-dominated Universe and a flat cosmological-constant- 
dominated Universe. These results will be presented in detail elsewhere . 

Figs. 1^ and ^estimate the probability of observing a set of parameters given an underlying 
model. Actual data will require solution of the inverse problem, estimating the probability 
of an underlying model given the data. We are currently exploring the accuracy with which 
all of the above parameters can be determined given a simulated data set ||l5| , pU| and to 
what extent parameter degeneracy in current experiments can be resolved . Preliminary 


results show that the true maximum of the likelihood function can be recovered with good 
accuracy from a parameter search routine, and that the errors in G approach the precision 
obtained here. 

The numerical results presented here demonstrate that G can be determined by realistic 
next-generation satellite experiments with a precision on the order of a few percent. Such 
a measurement will greatly solidify our knowledge of the gross properties of the Universe, 
have crucial bearing on the dark-matter and age problems, and will provide a stringent test 
of the inflationary hypothesis. 


6 







1 



0fwhm (degrees) 


FIG. 2. Same as Fig. 1, except for a reionized model with r = 0.5. 
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